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Aimlytical  methods  and  cooiputer  programs  have  been  developed  at  Douglas  for 
use  In  establishing  the  radiation  shielding  requirements  of  space  vehicle 
systems.  Current  techniques  for  performlrtg  space  radiation  shielding  analyses 
are  reviewed  in  the  context  of  the  digital  ccaaputer  programs  in  which  th^  have 
been  incorporated.  Calculations  dealing  with  the  transport  of  primary  charged 
particles  and  secondary  reaction  products  through  fflRtter  are  performed  by  the 
CHAKJE  progi^  for  geometrically  simple  systeme.  The  calculation  of  radiation 
levels  in  geometrically  complex  vehicle-shield  configurations  is  accomplished 
by  the  SICMA  program,  using  the  QUAD  generalised  quadratic  surface  geometry 
package.  Procedures  for  ascertaining  the  optimal  distribution  of  shield  mass 
subject  to  specific  dose  design  criteria  are  presently  under  development. 


1.0  IHTOODUCnOR 


Since  the  aaas  of  a  space  vehicle  shielding  systen  typically  represents  a 
considerable  payload  investnent,  an  accurate  sbnessaent  of  shielding  require* 
nents  is  essential.  The  overestlnation  of  shielding  requlreaents  can  lead  to 
an  excessive  vei^it  burden,  resulting  in  reduced  vehicle  performance  aarglns, 
increased  program  costs,  and  a  delay  in  space  exploration.  Conversely,  the 
underestinatlon  of  shielding  requlrenents  cam  seriously  Jeopaurdlze  alssion 
success . 

Studies  of  space  vehicle  systems  have  shovn  that  the  establishment  of  a  shield 
design  effectively  meeting  the  vehicle  and  mission  requirements  is  a  complex 
process.  Most  of  the  significant  factors  in  apace  vehicle  shielding  analyses, 
l.e.,  the  extamt  space  radiation  environment,  the  radiation  level  design 
criteria,  the  shield  system  gecanetric  framework,  etc.,  lure  sensitive  functions 
of  vehicle  and  mission  characteristics.  The  nature  of  these  relationahlpa  it 
Indicated  in  the  following: 

1)  Tte  mission  profile  governs  the  relative  importance  of  naturally 
occurring  earth  radiation  belts,  artificial  belts,  solar  radiation, 
and  galactic  cosmic  radiation;  the  dominant  Influences  are  due  to 
the  vehicle  trajectory,  the  mission  duration,  and  the  time  period 
of  the  mission  with  respect  to  the  solar  cycle. 

2)  Additional  influences  aw  due  to  those  ndsaion  requirements  affecting 
the  size  and  shape  of  the  crew  compartment  within  which  tolerable 
radiation  levels  are  to  be  maintained,  and  the  bloaedlcal  dose 
tolerances  established  for  crew  members 

3)  The  vehicle  configuration  is  Inportemt  since  the  internal  arraunge- 
mei’.t  of  work  and  rest  areas  controls  the  time-dependent  position  of 
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crev  nenbera;  the  location  of  fixed  equlpnent,  structure,  stores, 
etc.,  with  respect  to  work  stations  and  rest  6u«as  governs  their 
effectiveness  as  shielding  naterials. 

4)  Both  mission  and  vehicle  characteristics  influence  the  feasibility 
of,  and  benefits  from,  optimization  of  the  shield  design  by  dlstri> 
butlng  shield  mass  among  such  locations  as  crew  compartment  wall 
structure,  blowells,  localized  shields  and  portable  personal  shields. 

As  a  designer  and  manufacturer  of  space  vehicle  systems,  Douglas  has  actively 
punued  the  development  of  the  shielding  technology  required  to  integrate 
vehicle  and  mission  characteristics  into  shield  system  design.  The  shielding 
methodology  thus  developed  has  largely  been  Incorporated  in  digi.tal  computer 
programs  to  facilitate  its  application  to  radiation  protection  design  problems 
of  realistic  space  vehicle  systems.  These  design  methods  are  reviewed  here. 

Qnphasis  is  placed  on  the  shielding  technology  related  to  the  space  radiation 
environment,  rather  than  the  neutron  and  gaitsna  environments  from  on-board  nuclear 
power  systems.  Douglas  space  radiation  shielding  technology  is  discussed  in  two 
major  areas:  shielding  physics  and  applied  shielding.  The  shielding  j^yslcs 
efforts  have  concentrated  on  the  mathematical  formulation  of  the  transport  of 
radiation  through  matter,  usually  for  simple  geometric  systems.  Applied  shield¬ 
ing  has  desdt  with  the  calculation  of  radiation  levels  in  specific  source-vehlcle- 
shield  configurations,  with  the  associated  problems  in  treatrwnt  of  geometric 
parameters  and  shield  system  optimization. 
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2.0  SPACag  RADlAnOW  SHIEUDIWG  AJIALEIS 


Figure  1  llluetmtee  the  Interrelationship  of  several  cosiputer  programs  currently 
utilized  for  space  radiation  shielding  analyses.  The  manner  in  vhlch  a  shield¬ 
ing  optimisation  program  now  under  developwnt  will  be  integrated  with  existing 
programs  is  mxao  indicated. 

Topically,  the  radiation  analyses  requisite  to  space  vehicle  design  begin  with 
the  definition  of  the  extant  space  environment  in  space,  energy,  and  time.  Ilhe 
space  radiation  sources  which  are  impoirtant  to  space  syrtems  are  magnetically 
trapped  radiation,  solar  flares,  and  galactic  cosmic  rays.  A  computer  program, 
OGRE,  which  Douglas  has  developed  for  describing  geomagnet Ically  trapped  radia¬ 
tion,  presently  incorporates  the  latest  available  data  on  each  of  these  sources 
and  is  capable  of  determining  the  space  and  time  dependent  and  mission-integrated 
radiation  which  a  space  vehicle  will  encounter  on  a  specified  trajectory.  A 
complete  description  of  the  data  and  analytlca.l  techniques  used  In  the  OGRE  pro¬ 
gram  Is  given  in  Reference  1. 

Once  the  radiation  environment  has  been  defined,  computations  to  establish 
radiation  levels  can  be  performed,  the  calculatlonal  problems  becoming  increas¬ 
ingly  difficult  as  the  geometry  of  the  vehicle-shield  system  becomes  more 
complex  and  the  radiation  fields  become  more  extreme.  The  Douglas -developed 
CHARGE  program  is  the  primary  calculatlonal  tool  for  ascertednlng  beislc  dose 
attenuat''on  in  simple  geometries.  CHARGE  output  data  identifies  critical  organ 
dose  due  to  each  radiation  source,  l.e.,  trapped  protons,  trapped  electrons, 
solar  flare  protons,  secondary  electron-bremsstrahlurit:,  and  secondaury  nucleons 
from  proton  interactions  with  the  system  materials. 
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IWrE»HEy®IOSSHIF  OF  SPACE  RAniATIO!}  SHJELEm  PROGRWC 


PrCgTiUB 


Program  Function 


Mission  Specification; 

Pra>:tory 
Duration 
Laiinch  Date 


Specification  of  Shield 
Materials 


System  Description  j 

Vehicle  Configuratior 
Vehicle  Contents 
Man  Model  Geometry 
Time-deper.dent  location. 

of  Dose  Points 
Shield  Configuration 
Integration  Parameters 
for  Dose  Evaluation 


OGPE  HiOGRAM 

De-f^rmlnes  Mission  Integrated 
Radiation  Environment  Defined 
in  Intensity,  Energy,  aj.d  Time 


CHARGE  PROGRAM 

Computes  Dose  Transmission 
Throu^  Materials  in  a  Spher- 
icailly  Synsnetric  System 


SIGMA  PROGRAM 

Computes  Time-intecrated  Dose 
as  a  Function  of  Shield  Thick¬ 
ness,  Parameterized  Accordin' 
VO  Radi at in-  Type 


Printout  or 
Automat ic  Plot 


(UNDER  DEVELOPMENT) 


Optimization  Parameters  : 

Pote’^tlal  Shield  liocatiors 
Dose  Desigrs  Criteria 
Corstral'ts  Shield  Mass 
Co\cer.t at  Each 
hocat 

c’o  Criteria 


- - - 

SIGWA  Data  Parameterized 
Accordinp:  to  Each  Ray  Traced 
in  the  Numerical  Inteprrition 
Process 


r~z 


OPTIMISATION  PROGRAM 

Computes  Minimum  r^iield  Weight 

and  Optimal  Shield  Mass  Dis- 

tri  t.utio' 
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FIGURE  1 


Hie  Douglaa-'developed  SI9fA.  prcgrss  is  used  to  perform  dose  calcul&ti<  ns  for 
realistic,  geoisetrlcally  canplex  Space  vehicle  systems-  Hils  program  incor¬ 
porates  a  generalised  quadratic  surface  ^c»setry  package  by  means  of  which 
the  space  vehicle  is  defined  and  ray  tracing  operations  are  performed.  SICH4A 
computes  the  bioXogieal  dose  at  specified  receptor  locations  by  performing  a 
numerical  integration  of  angularly  dependent  dose  attenuation  kernels  over 
solid  angle.  Hie  dose  attenuation  keimels  necessary  for  this  arislysis  are 
customarily  constructed  by  curve-fitting  basic  CKABG2  dose  data 5  they  are 
then  evaluated  within  SICWA,  for  the  contputed  material  distribution  along 
each  ray  traced  in  the  angular  Integration  process,  to  obtain  the  dose  contri¬ 
bution  for  each  element  of  solid  angle. 

As  an  option,  the  dose  contribution  associated  with  each  ray  can  be  determined 
by  performing  s  complete  CHAB3E  radiation  transport  calculation  for  that 
particular  material  distribution.  Hiis  capability  la  illustrated  schematically 
in  Figure  1  by  the  dotted  arrcsw  between  the  SIGMA  and  CHARGE  programs  blocks. 

Hie  existing  OGRE,  CHARGE,  and  SIGMA  programs  will  be  supplemented  by  a  shield 
optimisation  program  cijrrently  under  development.  This  program  will  permit 
calculation  of  the  miRiaruia  shield  wei^t  meeting  specified  dose  design  criteria. 
For  this  purpose,  the  parametric  dependence  of  dose  or.  the  amount  and  dis- 
poBition  of  the  shield  mterlals  will  be  established  by  computations  largely 
perfonaed  with  existing  progreusa.  These  paraidetric  data  will  then  be  operated 
or.  by  an  iterative  technique,  based  or  the  Leigrsr.ge  method  of  indetermined 
multipliers,  to  deteniilne  the  opliaai  distribution  of  shield  zaaterisi  among 
several  candidate  locations  Ide.ntified  as  prograxp  iripu+y 


3.0  sHmpnc  PHTsics 


A  review  of  the  CHABGE  aethod  of  ccaputlng  the  transport  of  space  radiations 
through  aatter  is  given  In  this  section,  together  with  cosiparisons  of  CHARGE 
results  with  experlnental  aeasuresents  and  with  data  from  Independent  analyses.. 
Related  work  aimed  at  Improving  the  mathematical  models  and  basic  nuclear  data 
now  incorporated  In  the  CHARGE  program  is  also  presented. 

3.1  The  (mRGE  Program 

Shielding  j<x^3gram  CHARGE  (Reference  2)  provides  a  stradght-ahead  approximation 
solution  to  the  radiation  penetration  problem  for  primary  electrons  with  their 
associated  tremsstrahlung,  and  primary  protons  with  their  associated  secondary 
protons  and  neutrons,  for  geometrically  simple  shields.  Proton  and  electron 
analyses  use  a  modification  of  the  range-energy  relations  applied  to  as  many  as 
one  hundred  energy  groups.  Secondary  electron-brensstrahlung  transmlsaion  Is 
based  on  exponential  attenuation,  using  energy  dependent  non-elastic  cross 
sections  above  20  Mev  and  energy  independent  "removal"  cross  sections  below 
20  Mev. 


CHAJCJE  computes  the  dose  rates  and  energy  fluxes  as  a  function  of  penetration 
into  the  shield.  The  calculation  may  optionally  be  halted  before  the  entire 
shield  is  traversed  by  specifying  a  minimum  total  dose  rate  of  interest,  in 
which  case  it  proceeds  through  a  number  cf  shield  laml.na  Just  sufficient  to 
reduce  the  dose  rate  below  this  value. 


3.1.1  Primary  Particle  Ca-lcalatlons 


Thie  attenuation  of  prinsary  cn«riir»d  particles  is 
of  the  shield  materials.  The  primary  radiatlc." 


governed  primarily  by  ionisatlo 
dose  rates.  D  ,  are  y 

-  y  f 


nuioerical  evaluation  of  the  Integral 


~  E'(E,x)dE 


vbere  E'(E,x)  la  the  degraded  energy  at  shield  depth  x  (cm)  corresponding  to 
the  incident  energy  E,  and  is  calculated  from  tb*  range>energy 
relation,  (Mev) 

C  is  the  flux  to  dose  rate  conversion  factor,  (rad/hr)/ (Mev/cm^-sec) 

^(E)  is  the  incident  inward  directed  primary  particle  flux,  (paz^lcles/Mev- 
cm  -sec) 

/„^(E')  is  the  macroscopic  cross  section  for  non-elastic  collisions, 

(«‘^) 

E^  is  the  maximum  energy  of  the  incident  particles,  (Mev) 

E^(x)  is  the  energy  >rtiich  has  range  x  in  the  shield  material,  (Mev;. 


llie  incident  spectrum  is  divided  into  energy  groups  and  the  'shield  materials 
into  layers  for  the  penetration  analysis.  In  contrast  to  usual  multi-group 
treatments,  the  spectrum  shape  is  retained  within  each  energy  group  when 
performing  the  dose  integral.  Ihls  approach  yields  much  more  accurate  results 
than  is  possible  using  many  tiroes  the  number  of  energy  groups  with  a  constant 
spectrum.  Since  the  mean  free  path  for  non-elastic  coliisions  is  large  com¬ 
pared  to  rsual  shield  thicknesses,  tnd  varies  little  over  wide  ranges  of  energy, 
the  integral  in  the  nuclear  attenuation  term  is  closely  approximated  by  a 
suirmation  using  average  group  energies.  Hie  incident  group  of  jharged  partici  .s 
Is  attenuated  through  the  shield  by  an  integration  of  ^  E'(E,x),  determining  the 
group  boundary  energies  as  a  function  of  x.  The  treatment  of  electrons  and  protons 
differs  slightly  because  of  their  different  range-energy  c'  iracteri sties. 


Charged  fttrtlcle  Spectra  -  llhe  ambient  electron  and/or  proton  spectra  may  be 
specified  eis  differential  spectra  in  tabular  or  functional  form.  '3ie  energy 
dependence  of  these  spectra  is  almoi’t  unlimited,  and  the  energy  ranges  con- 
sldered  are  limiter'  only  by  the  computer  capabilities.  Uiese  Incident  charged 
particle  spectra,  which  may  be  viewed  as  essentially  continuous  functions  of 
particle  kinetic  energy,  are  divided  into  energy  bands  or  groups,  the  number 
(<  lOO)  and  spacing  of  irtiich  are  controlled  by  input  data.  Each  of  these  groups 
may  be  completely  described  by  its  upper  and  lower  energy  limits  and  the  energy 
variation  of  the  spectrum  within  the  group,  with  CHARGE  ensuring  that  the 
spectrum  within  any  ^roup  is  a  single,  continuous  function  of  the  energy- 

Functional  forms  allowed  for  the  ambient  differential  spectra  Include  (l)  a  power 
law  or  (2)  an  exponential  dependence;  the  spectrum  may  be  specified  for  each  of 
several  energy  ranges.  Input  data  for  either  analytic  form  may  consist  of: 

1)  Total  ambient  particle  densities. 

2)  Power  of  E  or  multiplier  of  E  in  exponent  over  each  energy  range 
(used  for  interpolation  to  energies  not  given  in  input). 

3)  Vedues  of  differential  spectrum  at  input -sped fled  energies. 

U)  '  onstants  which  multiply  the  power  or  exponent  of  E  to  give  the 
density  at  each  energy. 

Proton  Energy  Loss  Rates  -  The  proton  range-eriergy  computation  is  based  on  the 
relativistic  energy  lose  rate  for  "heavy”  particles  traversing  a  material  and 
losing  energy  by  10’;''zation  (l.e  ,  elastic  collision  with  atomic  electrons). 

This  formula,  give”^-  in.  many  references  (e.g.,,  Reference  3)>  Is 


where  e  is  the  electronic  charge, 
m  Is  the  electronic  raass, 

Z  is  the  charge  of  the  incident  particle,  units  of  e 
V  is  the  velocity  of  the  incident  particle, 

N  is  the  density  cf  stopping  atoms,  and 
B  is  the  "stopping  number". 


According  to  the  Bethe-Bloch  theory 


B  =  Z 


fci(l-/8^)  -  6 


where  Z  is  the  atomic  number  of  the  stopping  material,  I  is  the  average  ioniza¬ 
tion  potential  (energy  to  remove  one  electron)  for  the  stopping  atoms,  is  \ 
divided  by  c,  the  velocity  of  light,  and  6  is  the  “density  effect”.  The  density 
effect,  6,  is  evaluated  by  the  empirical  form  given  by  Sternheimer  in  References 
h  and  5* 


Proton  Range-Energy  Calculation  -  The  proton  range  is  defined  as  the  distance 
penetrated  before  coming  to  rest.  For  a  proton  of  energy  incident  on  a 
shield  layer  of  thickness  x,  CHARGE  computes  the  transmitted  energy  E^  and  the 
range  corresponding  to  the  transmitted  proton  ener^  R(EL): 

c: 


P 


c 


i/y 


R 

0 


y 


)0 


To  accoiEpli«*h  this,  ar*  initial  eatiaate  of  the  transmitted  energy  is  made 


E,  =  E  +  X  -s-i 
1  o  dx 


and  the  approximate  energy  loss  rate  dE^^/dx  established.  Assuming  the  range- 
energy  relationship  to  be  of  the  form 


dE  dE  o 
dx  *  ^  Tr 


Uie  values  of  7  and  required  to  establish  and  RCE^)  are  obtained  by 

substituting  E^  and  E^^  in  the  latter  expression  emd  rearranging  terms: 


!2  /‘ 

■  y  /  < 


y  =  1  + 


|n(-dE^/dx)  -  Ih(-dEydx) 


Electron  Energy  Loss  Rates  -  CHARGE  uses  the  average  energy  loss  formula  for 
relativistic  electrons,  including  exchange,  spin,  density,  and  quantum  mechani¬ 
cal  effects,  and  an  approximate  correction  for  energy  loss  due  to  bremsstrphlung 
emission: 


dE  Une 


vbere 


This  fomula  results  from  combining  several  eijuatlons  given  by  Bethe  and  Asbkln 
in  Reference  3*  At  energies  below  a  specified  value  E^,  dE/dx  is  held  constant, 
as  for  protons. 


Bremsstrahlung  production  is  of  little  importance  as  an  energy  loss  laechanisin 
for  electron  energies  normally  encountered,  edthough  it  is  a  significant  source 
for  dose  determinations. 


Electron  Transmission  Calculation  -  Electron  transmission  is  computed,  option¬ 
ally,  by  a  rauige-energy  method  or  by  semi -analytical  transmission  factors.  Hie 
electron  range-energy  option  assumes  the  same  functional  relationship  between 
penetration  and  energy  as  is  used  for  protons.  The  semi-analytical  transmission 
factors  are  derived  from  Monte  Carlo  data  (Reference  6);  the  transmission  factor 
treatment  is  normally  used  in  preference  to  the  reuige-energy  method  for  shields 
of  a  single  material.  In  order  to  improve  results  obtained  through  the  use  of 
transmission  factors  for  laminates  of  various  materials,  electron  penetration 
studies  using  an  analytic  trauisport  code  are  in  progress.  The  data  from  these 
investigations  will  be  incorporated  in  CHARGE  ai'  they  become  available. 


3.1.2  Secondary  Radiation  Calculations 


Secondary  radiations  treated  by  CHARGE  are  electron  bremsstrahlung  and  neutrons 
and  pn  tof.s  from  primary  proton  interactions  with  shield  materials..  The  trans- 
piTt  fu:  ctlor.s  for  secoridary  protons  are  identical  to  those  for  primary  protons 


f  the  same  er.ergy 


a©  secondary  particle  dose 


rates  are  represented  by  the  Integral: 


la  the 


T/»U-x'.y 


cross  section  for  production  of  secondary 
particle  type  fi  with  energy  from  prlaary  type  a 
irtth  the  energy  (cii“^) 

Is  the  transport  function  for  particle  type  fi,  created 

with  energy  ^  «t  position  x'  to  the  detector  at  x, 
(dimensionless) 

))  li  the  flm  to  doee  rate  eooverelon  for  ptttlcle  tjpe  fi 
«t  deeradod  energr  E'o,  (nS/lSL 

\» 


iHsv/ca^-aec, 


— -  Oroe,  .action  det.  for  prodoctlon  of  „cond«ry  nucl«,„.  hf 
nucleon.  .trlkln«  cocple,  nuclei  b.,e  been  generated  end/or  compiled  et  DougUe, 
-  reported  In  Bef.renc,  T.  IM.  tre.t«„t  eeperete.  the  pnniuct  nucleon,  into’ 

.  fon^rd.peehed  hlgh-enerar  (cMc.de)  ph...  end  en  Uotioplc  loo-energy  (.«,or.. 
•-ion)  phMe.  CKAICE  Incorporete.  the  ,.cond«y  nucleon  production  dete,  .. 
pr..ent.d  In  Reference  7,  «th  the  .ingle  exception  thet  11  cede  pertlcle. 

treated  m  being  .itted  etrlght-ehe^l;  AlMlller  et  1  he,.  «.<*,  thl, 
—UMptlon  to  be  .lid  for  u  .,c„ng.ay  nucleon,  (Reference  8).  i„  eSAX^, 

..conder,  nucleon.  Me  clcul.ted  for  .  .i„*ie  cherecterl.tlc  .l.«nt  In  ealh 
shield  material. 


Since  the  low  energy  protons  are  much  more  effective 
are  high  energj^  protons,  this  source  la  treated  mere 


in  depositing  energy  than 
accurately.  For  low  energy 


secondary  protons,  the  source  is  computed  by  numerically  integrating  the  pro¬ 
duction  cross  section  over  each  shield  region  thickness.  Ihe  energy  dependence 
of  the  priEHiry  protons  is  treated  by  assigning  an  average  energy  E^(x’)  over 
each  energy  group  i.  However  the  variation  of  Ej^(x')  with  penetration  depth  x* 
is  accurately  cooputedi  so  taat  region  thicknesses  x  are  not  restricted  to  small 
values.  The  secondary  source  energy  ^(x'}  is  also  treated  as  position-dependent. 

For  high  energy  secondary  protons,  the  dependence  of  both  primary  and  secondary 
proton  ener<ri**s,  E^(x')  and  Ej^(x'),  on  position  within  the  region  is  not  treated. 
They  are  considered  to  be  constant  over  the  region  thickness  x,  because  of  their 
large  range. 

The  treatment  of  the  spatial  distribution  of  secondary  neutron  production  is 
similar  to  that  used  for  high  energy  secondary  protons.  Computed  values  of 
total  doee  are  insensitive  to  the  detailed  treatment  of  neutron  production 
because  of  the  significantly  lesser  attenuation  of  neutrons  along  the  trans¬ 
mission  path  as  comptired  to  protons. 

Secondary  neutrons  are  attenuated  using  an  energy-dependent  (group-wise)  removal 
cross  section  for  neutron  energies  greater  than  20  Mev.  Bie  neutron  non-elastic 
cross  section  data  for  this  computation  aure  presented  in  Reference  7*  For  lower 
neutron  energies,  a  constant  valtiC  of  the  removal  cross  section,  based  on  the 
shield  material  traversed,  is  applied. 

Electron  Ere  jasst rah lung  -  Several  formulations  of  the  bremsstrahlung  production 
cross  section  are  Incorporated  in  CHARGE,  providing  a  meauns  of  comparing  results 
fn.iK  the  various  techniques  for  calculating  the  source  of  this  dcse  component.; 
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Bie  eoplrlcal  treatment  by  Evana  (Reference  9)  i*  the  aimpleat  model  and  ylelda 
reaulta  idilch  agree  reaaonably  aell  with  a  combination  of  Bcthe-Heitler  non- 
acreened  crosa  aectlona  and  varloua  croaa  aectlona  compiled  by  koch  and  Mot* 
(Reference  10).  The  photona  thua  generated  are  tranaported  from  the  aource 
point  to  the  detector  location  by  exponential  attenuation  and  buildup.  The 
calculated  photon  apectra,  which  are  tranaadtted  throu^  the  ahleld,  bracket 
the  experimental  apectra.  of  Edelaack  (Reference  11 )  for  1.0  to  2.0  Mev  electrona. 

An  experimental  program  haa  been  initiated  by  Douglaa  to  determine  enei’gy  and 
angular  dlatributlona  of  photona  produced  by  electrona  over  the  energies >  and 
in  the  mateirlala,  applicable  to  apace  vehicle  shielding. 

3.1.3  Dose  Rate  Calculations 

Charged  Particlea  -  At  each  ahleld  region  interface,  CHARGE  cooputea  the  dose 
rate  due  to  incident  charged  particlea,  including  the  effects  of  absorption 
and  energy  loss  in  the  region  already  traversed.  Since  the  doses  are  not 
computed  at  each  mesh  interval,  the  nuniber  of  regions  is  determined  by  the 
desired  number  of  dose  points  and  the  actual  number  of  phyalcally  distinct 
layers,  while  the  number  of  mesh  intervals  is  determined  by  error  limits  on 
the  accuracy  of  the  range-energy  integration  and  ccmputer  running  time  iimlta- 
tlona.  In  general,  the  mesh  intervals  should  be  thinner  than  ten  per  cent  of  the 
range  of  the  highest  energy  pexticles. 

Vo  evaluate  charged  particle  contributions  to  biological  dose,  CHARGE  currently 
has  two  options:  one  of  these  assumes  the  "relative  biological  effectiveness" 
(RBE)  of  deposited  energy  is  a  function  of  particle  energy,  the  other  assumes 
RBE  is  a  function  of  the  energy  loss  rate  or  linear  energy  transfer  (LET). 


3.2.1  Breasstrahlung  Spectra  in  Alumlniai 

In  C3IARGE,  electron  brenitrahlung  production  may  be  calculated  by  the  Bom 
approximation  or  by  a  almple  formula  given  by  Kulenkampf  (Reference  12). 

Figure  2  preaenti  a  compcurlion  of  bremsitrahlung  photon  spectra  calculated 
by  each  of  these  methods  vlth  experimental  data  of  Edelsack  et  al  (Refer¬ 
ence  U).  These  data  have  been  renormallted  to  represent  the  photon  flux 
Integrated  over  all  solid  angles.  The  spectra  »Te  those  transmitted  thmugh 
an  aluminum  slab  Just  thick  enough  to  stop  the  monoenergetlc  Incident  electrons.; 
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FIGURE  3 


DOSE  CALCUUTiON  FOR  ARTIFICIAL  ELECTRON  BELT 


FIGURE  3 


DOSE  IN  WATER  FROM  A  TYPICAL  VAN  ALLEN 
PROTON  SPECTRUM  TRAVERSING  A  SPHERICAL  ALUMINUM  SHIELD 


since  the  Kulenkaapi'  (or  "K")  nodel  predicts  identical  relative  spectra  for  all 
electron  energies,  aultlple  curves  appear  only  where  the  photon  absorption  differs. 
The  K  option  is  normally  used  for  shielding  analyses  because  It  predicts  a  spectrum 
in  better  jgreement  with  experimental  data  for  a  wide  range  of  materials,  and 
because  it  is  a  very  simple  form. 

3-2.2  Dose  Calculation  for  Artificial  Electron  Belt 

Dose  cadculations  were  performed  for  an  electron  spectrum,  ai>praxisntely  repre¬ 
senting  the  artificial  electron  belt,  incident  on  an  infinite  aluminum  shield. 

Bie  electron  transmission  calculation  used  the  range-energy  option  (52  electron 
groups)  rather  than  the  more  accurate  semi-analytic«LL  transmission  factor  sethod. 
The  bremsstrahlung  source  was  computed  from  empirically  corrected  Bom  approxl- 
nation  cross  sections. 

Results  are  compared  with  data  published  by  Russak  (Reference  13)  in  Figure  3- 
The  agreement  is  generally  good.  Differences  in  the  two  electron  dose  curves  may 
be  attributable  tc  the  use  of  different  flux-to-dose  conversion  fact'^rs. 

Ihese  data  are  indicative  of  the  dose  rate  which  would  be  received  by  an  astronaut 

performing  operations  external  to  the  space  vehicle  and  therefore  protected  only 

by  a  space  suit.  The  dose  rate  under  these  conditions  is  relatlvel.-  insensitive 

dE 

to  the  Incident  electron  energy  spectrum,  since  ^  is  not  a  strong  function  of 
energy  over  the  electron  energy  range  of  interest., 

3.2.3  Proton  Dose  Comparison  with  Analytic  Range-Energy  Data 

Corry  and  Stogryn  ctdculated  the  dose  behind  an  aluminum  slab  exposed  to  a 
monodlrectlonal  proton  flat  typical  of  the  Van  Aller.  belt  (Reference  l4). 


The  analytic  calculation  assuioed  the  proton  range  to  be  approximated  by  a  single 
power  law  of  energy  for  all  protons.  CHABGE  results,  obtained  using  the  properly 
adjusted  isotropic  ambient  flux,  are  compared  with  this  work  in  Figure  4.  T^ie 
agreement  is  excellent. 

The  cusp  cf  the  cvirve  is  a  direct  result  of  the  fact  that  the  incident  spectrtm 
was  specified  only  for  energies  above  10  Mev.  As  the  lowest  energy  protons  slow 
down,  their  biological  effectiveness  increases,  this  effect  predominating  until 
the  penetration  corresponding  to  the  range  of  10  Mev  protons  is  reached  (cusp  of 
figure),  at  which  point  proton  absorption  predominates  and  dose  decreases  with 
depth.  This  tjiJe  of  behavior  appears  whenever  the  incident  spectrum  is  unspecified 
for  low  energies. 

j.2.U  Dose  Calculation  for  Van  Allen  Protons 

A  typical  Van  Allen  proton  spectrum  was  represented  for  two  energy  ranges  by 

(30<E<80>fev) 

^  (8o<E<700  >tev) 

and  normalized  to  a  flux  of  one  proton/ cm  -sec.  The  detailed  dose  conirtbutlons 
computed  by  CHARGE  for  primary  protons  and  secondary  protons  and  neutrons  behind 
an  aluminum  shield  are  shown  in  Figure  5-  A  ccnparlson  is  made  with  results 
from  an  ionization  calculation  which  assumes  that  all  energy  losses  are  due  to 
ionization  since  the  attenuation  of  primary  charged  particles  is  governed  prin¬ 
cipally  by  this  method.  For  this  pi^Dton  spectrum,  the  ionization  method,  which 
requires  much  less  computer  time,  yields  a  good  approximation  to  the  total  dose 
for  small  shield  thicknesses,  but  is  conservative  with  increasing  shield  depth. 


3'2.5  Doge  Cadculatlon  for  Solar  Pliire  Protoim 


In  Figure  6,  the  doge  received  from  a  golar  flare  of  the  November  12,  19^0  type 
ig  pregented  as  a  function  of  aluminum  ahield  thickneaa.  The  differential  proton 
gpectrum  vaa  approximated  in  three  energy  regiona  by 

^  -  1.77  X  10^3  E‘3  (30<E<80Mev) 

^  =  9-62  X  E"^  (80<E<4ii0  Mev) 

^  =  6.63  X  10^®  E"^’^  (W^O<E<6600  Mev) 

irtiich  coiresponds  to  an  integrated  flux  of  9-0^  x  10^  protor.s/os'^ .  The  Ioniza¬ 

tion  cadculatlon  gives  doses  within  ten  per  cent  of  the  detailed  calculation  for 
shield  thicknesses  less  than  30  g/cc  *  For  shield  thicknesses  greater  than  about 
30  g/cm  ,  the  dose  from  secondary  radiation  is  important  because  of  the  rapid 
attenuation  of  the  primary  proton  dose.  This  behavior  is  characteristic  of  the 
steep  primary  proton  spectra  typically  associated  with  solar  flare  evt  its. 
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DOSE  PER  FLARE  (RAO) 


Dost  IN  ^'MER  FROM  NOVEMBER.  I%0 
SOLAR  FLARE  PROTONS  TRAVERSING  A  SPHERICAL 

ALUMINUM  SHIELD 
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FIGURE  6 


4.0  APPLIED  SHIELDING 


The  evaduation  of  shielding  requirements  for  a  space  vehicle  system  necessitates 
a  detailed  treatment  of  the  actual  system  geometry.  Dose  computations  are  orten 
extremely  sensitive  to  inhonogei  cities  in  the  spatial  dib‘ributlon  of  structural 
and  shield  materials,  eq’uipment  and  stores,  etc.,  so  that  an  aeoirate  description 
of  the  material  distribution  is  essential.  Since  the  shielding  physics  calcula¬ 
tions  of  radiation  -  materials  interactions  are  restricted  to  simple  geometries, 
these  methods  are  normally  applied  to  construct  kernel  functions  which  can  then 
be  used  for  dose  evaluations  by  applying  ray  tracing  procedures  to  complex  geo¬ 
metric  systems.  At  Douglas,  the  geometric  description  of  a  realistic  space  vehicle 
is  performed  by  the  SICMA  program,  which  incorporates  a  numerical  procedure  for 
angular  integration  of  dose  attenuation  kernels  constructed  by  the  CHARGE  program. 
Geometric  operations  within  SIGMA  are  performed  using  the  C^AD  generalized  quadric 
surface  geometry  package. 

It  can  be  Important  to  consider  carefully  the  geometric  distribution  of  biological 
materials  about  the  critical,  radiation-sensitive  body  orga;  s  of  crew  members. 

For  thinly  shielded  systems,  particularly  for  extravehicular  operations  perfoimied 
in  a  space  suit,  the  shielding  afforded  by  body  tissues  may  represent  the  greater 
part  of  the  total  effective  shielding.  Therefore  studies  comparing  dose  computa¬ 
tions  for  various  extant  man  models  have  been  performed.  11113  work  has  resulted 
in  an  improved  mathematical  model  being  outlined,  which  is  cased  on  a  Cartesian 
coordinate  representation  of  body  material  distributions  Hie  objective  of  this 
work  is  +0  provide  a  stsindard  against  which  results  from  simpler  models  can  be 
measures. 

The  detenui nation  of  shield  weight  requirements  for  specified  vehicles  and  missions 
is  str'^ngly  dppprdt"^t  t  the  allowable  shield  gpo-oetric  framework.:  Hie  dlstributio 


of  shield  materials  among  various  locations  onboard  the  vehicle,  and  the  shaping 
of  shield  thicknesses  over  extended  areas,  can  be  Isipoirtant  vhen  seeking  to  meet 
the  radiation  level  design  criteria  for  alniinian  shield  system  wei^t.  At  present, 
Douglas  is  adapting  paverful  mathematical  optimization  techniques,  previously 
used  in  neutron  and  ganna  shielding  studies,  for  ajjpli cation  to  space  radiation 
shielding  optimization  investigations. 

4.1  Detailed  Treatment  of  Vehicle  Geometry 

The  treatment  of  actual  vehicle-shield  geometries  in  radiation  level  calculations 
involves  the  use  of  ray-tracing  techniques.  The  geometric  description  of  physical 
regions,  and  the  determination  of  path  lengths  traversed  by  radiations  penetrating 
the  vehicle  regions,  is  accomplished  by  the  QUAD  geometry  package.  The  QUAD  sub¬ 
routine  is  under  the  control  of  the  SIGMA  program  which  sets  up  the  angular  inter¬ 
vals  for  a  numerical  integration  of  dose  contributions  over  the  solid  angle  about 
each  receptor  locations,  and  i^ich  controls  variations  in  shield  thicknesses, 
shield  materials,  etc.,  for  parametric  shielding  studies. 

4.1.1  Geometry  Calculations 

The  Douglas  QUAD  subroutine  is  a  powerful  and  efficient  method  for  describing 
geonetric  systems  by  generalized  quadric  surfaces.  The  system  geometry  is 
described  by  .i  set  of  material  regions  with  one  or  more  boundaries  which  are 
defined  by  equations  of  the  form: 

Aq  +  A^x  +  A^y  +  A^z  +  Aj^x^  +  A^y^  +  A^z^  +  A.^y  +  Agyz  +  A^xz  =  0 

where  ...»  are  input  constants  and  (x,  y,  z)  are  the  usual  Certeslan 

coordinates.  The  version  of  this  subroutine  used  with  SIGMA  allows  up  to  100 
surfaces  and  100  regions  for  system  description.. 


cyjAD  performs  the  following  functions  under  the  control  of  the  SIGMA  program: 

1)  Input  of  the  system  geometry  description, 

2)  Determination  of  the  region  occupied  by  a  point  with  specified 
coordinates. 

3)  Determination  of  the  regions  traversed,  the  path  length  in  each 
region,  and  the  boundary  crossed  in  leaving  each  region,  for  a 
ray  specified  as  to  origin  and  direction. 

QUAD  incorporates  several  features  which  ease  input  requirements,  reduce  the 
possibility  of  input  and  calculational  eiTors,  and  promote  computational 
efficiency.  Hiese  include  the  follovring: 

1)  The  ’’ambiguity  index”  is  internally  computed  for  each  boundary  of 
each  region.  This  index  has  the  value  +1  or  -1  depe.^ding  on  the 
location  of  the  region,  inside  or  outside  respectively,  with  respect 
to  the  surface  forming  this  boundary. 

2)  Surfaces  are  defined  separately  from  the  geometric  regio^is  to  eliminate 
duplication  of  data  when  a  surface  forms  a  boundary  of  more  than  one 
region . 

3)  All  dislarices  to  s’orfaces  are  computed  and  stored  independent  of  the 
region  ooundaries,  thus  elimiriating  repetitive  calculation  of  such 
quantities  for  those  surfaces  which  bound  a  r.umber  of  regions  along 
a  specified  ray., 

.  i  .  I\)so  Kvaluation 

D.e  'I  ’M\  progrfir  calculates  tne  dose  received  at  each  of  a  set  of  specified 
p  ♦  :v  i  f'T  al  to  a  vehi clo-sMeld  system  by  a  r^.imerical  Gimpsor^s  mle 


integration  of  the  equatlo;,; 


K(7,^d0cl// 


where  r  =  dose  point  location 

D(r)  =  dose  received  at  point  r 
ft  =  the  cosine  of  the  polar  angle 
0  =  the  asimuthal  angle 

^  =  y  1-^  (cos  Oi  +  sin  )  (a  unit  vector  in  the  direction 

defined  by  H  and  O) 


The  numbei  of  intervala  in  the  U  ®  integrations  are  variable,  with  systems 
of  pronounced  asynmetry  requiring  e  larger  number  of  intervals.  The  kernel 
K(?,^  represents  the  dose  that  would  be  received  at  r  if  the  material  thickness 
encountered  along  a  ray  from  r  in  the  direction  of  H  were  spherically  symnetric 
about  r.  These  material  thicknesses  are  calculated  by  the  Douglas  compiter  sub¬ 
routine  QUAD;  the  kernel,  K(r,S)  is  evaluated  for  the  computed  material  distribu¬ 
tion  alc^iHI  the  ray  from  either  curve-fit  CHARGE  data  or  from  a  complete  CHARGE 
calculation 


TrfO  options  which  can  be  used  concuirently  are  provided  for  the  calculation  of 
the  dose  received  throur^h  various  outer  sections  of  the  vehicle.  The  first 
provides  a  dose  breakdowi’  according  to  specified  geometrical  regions  through 
which  the  radiation  enters  the  vehicle,  and  tne  second  limits  the  ranges  of 


integratior  of  the  angular  variables.: 


integration  of  the  equation: 


K(?,ddOd// 


where  r  =  dose  point  location 

D(r)  =  dose  received  at  point  r 
H  =  the  cosine  of  the  polar  angle 
0  =  the  asimuthal  angle 

(cos  Oi  +  sin  S,])  ^  Mk  (a  \tnlt  vector  in  the  direction 

defined  by  U  and  O) 


Hie  nunbei  of  intervals  in  the  n  and  0  integrations  are  variable,  with  systems 
of  pronounced  asymnetry  requiring  a  larger  number  of  intervals.  Hie  kernel 
K(?,^  represents  the  dose  that  vould  be  received  at  r  if  the  material  thickness 
encountered  along  a  ray  from  r  in  the  direction  of  i}  vere  spherically  symmetric 
about  r.  These  material  thicknesses  are  calculated  by  the  Douglas  compiter  sub¬ 
routine  QUAD;  the  kernel,  K(ry5)  is  evaluated  for  the  computed  material  distribu¬ 
tion  alcw  the  ray  from  either  curve-f'it  CHARGE  data  or  from  a  complete  CHARGE 
calculation. 


Tvo  options  which  can  oe  used  concun*ently  are  provided  for  the  calculation  of 
the  dose  received  through  various  outer  sections  of  the  vehicle.  The  first 
provides  a  dose  breakdown,  according  to  specified  geometrical  regions  through 
which  the  radiation  enters  the  vehicle,  and  the  second  limits  the  ranges  of 
integratior  of  the  angular  variables. 


Pr^ 

Asi( 
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Frognm  SIGNA  Incorporates  several  features  to  facilitate  panusetric  studies. 


Aaong  these  are  the  follovlng: 

1)  Multiple  orbits  are  treated  slaultaneously  during  the  geoastrle 
calculations.  Each  orbit  amy  entail  several  radiation  sources 
and  the  dose  received  frosi  each  prlmry  and  secondary  radiation 
cosiponent  is  printed  out. 

2)  Multiple  iBoterial  sets  are  alloved,  and  as  a  result,  each  aaiterlal 
region  can  have  several  different  densities;  this  allows  a  deter- 
ndnatiai  of  the  variable  shielding  effects  of  equipnent,  supplies, 
etc.,  and  facilitates  paranetrlc  studies  of  shield  thickness  where 
geometric  effects  are  negligible  and  changes  in  shield  thickness 
can  be  approximated  by  changes  in  material  density. 

3)  The  capability  to  input  multiple  material  sets  pexmits  dose 
attenuation  data  for  various  shield  materials  to  be  approximated 
by  expressing  material  thickness  in  terms  of  an  equivalent  density 
of  the  base  material.  'Rils  approximation  is  acceptable  for  proton, 
electron,  and  x-ray  or  gamom  ray  doses  when  the  materials  are  not 
very  dissimilar  and  the  attenuation  is  not  too  great.  For  large 
parametric  studies,  it  has  been  found  that  this  can  save  a  great 
deal  of  computer  time  without  appreciably  affecting  the  accuracy 
of  the  results. 

4)  Multiple  material  seta  and  multiple  orbit  sets  are  all  treated  during 
a  single  angular  integration,  so  that  geometric  computations  are  not 
repeated  unnecessarily. 

5)  Dose  transmission  may  be  calculated  by  several  alternate  methods 
including  a  basic  transport  calculation  for  the  material  distribution 
along  each  ray  traced.  For  example,  a  complete  CHARGE  calculation 


can  be  perforaed  tor  eadi  rey  traced.  These  calculations  can 
include  the  effects  of  anisotropic  source  distributions. 

6)  CHADCaS  data  for  kexml  evaluations  are  input  in  tabular  fora  and 

curve«fitted  internally  by  exponentials  or  power  lavs  over  each 
interval.  This  e<^iderably  siapUfies  the  preparation  of  input 
data  over  aethods  rsfuire  the  deteraination  of  coefficients 

for  the  kenteli  external  to  the  prograa.  Anisotropie  sources  can 
be  treated. 

7)  Computed  dose  data  tor  various  dose  receptor  locations  in  the 
vehicle  can  be  tiae-wel^ted  according  to  prospective  work-rest 
cycles  and  suuaed  to  provide  representative  total  dose  data  for 
an  entire  aission.  Such  data  are  obtainable  as  s  function  of 
shield  thickness  when  the  option  to  use  successive  sets  of  material 
data  for  each  shield  region  Is  exercised. 

8)  SICMA  output  data  on  computed  doce  as  a  function  of  such  parameters 
and  variables  as  shield  thickness >  orbital  altitude,  orbit  inclina¬ 
tion,  etc.,  can  be  plotted  directly  on  an  G-C  U020  high  speed  niero- 
fiJa  recorder  operating  as  peripheral  equipamnt  on  the  Douglas  IBH 
7094  costputer.  Considerable  flexibility  Ir  arrangement 

labeling  of  output  data  is  possible. 

1*.2  Related  Geoaietrlc  Considerations 

In  addition  tn  the  detailed  treatment  of  vehicle  geometry,  there  exist  other 
considerations  of  a  geometric  nature  vhlch  can  be  important  in  establishing 
shielding  requirements  for  a  given  space  vehicle  end  mission.  One  such  con¬ 
sideration  is  the  geometric  representation  of  astronauts,  particularly  the 
spatial  distribution  of  biological  avaterials  shielding  critical  body  organs. 
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Another  Is  the  eBtablishiient  of  the  gecaaetrlc  frwneworK  within  which  shield 
naterials  are  positioned;  l.e.,  the  determiiMition  of  an  optlaal  distribution 
of  shield  aass  is  subject  to  the  ruaber  and  type  of  variables  defining  the 
potential  geonetrlc  locations  of  the  attenuating  oaterials. 

4.2.1  Doeiaetrlc  Anatoadcal  Model  of  Man 

A  significant  factor  in  spaceoradiation  analysis  is  the  shielding  worth  of  body 
tissues  enclosing  the  radiation-sensitive  body  organs.  Thus,  one  of  the  iinport- 
ant  models  in  computer  programs  designated  for  analyzing  space-radiation  shield¬ 
ing  requirements  is  an  anatcmdcal  model  of  man.  Several  radiobiological  man 
models  of  varying  sophistication  have  been  reported  in  the  space-radiation 
literature.  Some  of  these  models  are  decided  over-simplifications  of  the  human 
anrtomy  and  of  critical  organ  locations;  other  models  treat  the  spatial  distribu¬ 
tion  of  body  materials  in  considerable  detail.  Most  of  the  extant  man  models 
are  deficient  in  combining  accuracy  in  the  representation  of  body  tissue  distri¬ 
butions,  for  a  variety  of  potential  body  positions,  with  efficiency  in  computer 
application  of  the  model  and  in  preparation  of  the  anatomical  input  data 
requisite  to  its  use. 

Comparison  of  computed  dose  results  for  two  of  the  comnonly  applied  man  models, 
positioned  within  a  realistic  space  vehicle,  are  shown  in  Figure  7  and  8.  Fixed 
vehicle  structure  and  equipmsit  are  worth  approximately  2.2  g/cm  .  The  "average 
tiiisue  depth"  model  assumed  0.3  cm  of  tissue  shielding  over  a  2n  solid  angle 
and  Infinite  shielding  over  the  remainder.  The  elliptical  cylinder  model  is 
comprised  of  two  elliptical  regions,  one  representing  the  head  and  neck,  and 
other  the  trunk:  cylinders  were  dimensioned  according  to  mean  dimensions  of 
USAF  flying  personnel  (Reference  15).  Calculations  based  on  this  model  are  for 
a  dose  poirit  0.,3  cm  below  the  surface  of  the  smaller  cylinder,  at  eye  level 
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Differences  betwen  computed  dose  results  for  the  two  sKidels  are  a  direct  result 
of  the  assumption  of  2«  Infinite  shielding  for  the  average  tissue  depth  model. 
For  the  readily  attenuated  trapped  electron  dose  contribution,  this  assumption 
results  in  an  over-estimation;  i.e.,  the  shielding  effect  is  act’ially  greater 
than  2n  because,  in  the  elliptlceO.  model,  the  slant  path  along  rays  deviating 
significantly  froci  the  normal  is  gx*eater  than  the  average  tissue  depth  of 
0.3  cm.  For  the  more  penetrating  protons  and  electron  bremsstrahlung,  the 
assumptions  of  2«  infirdte  shielding  results  in  an  under-estlmation  of  dose, 
ms  is  because  the  shielding  offered  by  the  head  is  not  effectively  infinite, 
as  it  is  for  electrons.  Comparable  results  were  obtained  for  dose  to  the  skin 
and  to  blood-forming  organs. 

In  the  comparison  shown,  the  elliptical  cylinder  model  was  taken  to  be  the 
standard;  results  based  on  this  model  are  considered  "correct."  However,  the 
actual  correct  values  are  not  known  Without  resort  to  either  empirical  data 
or  still  more  sophisticated  mathematical  models.  Shield  weight  requirements 
are  quite  sensitive  to  the  man  model  used,  as  indicated  by  the  pairs  cf  shield 
thicknesses  shown  for  each  of  the  two  design  dose  criteria;  uniformly  distributed 
shielding  weighs  3500  lb  per  g/cm  for  the  reference  configuration.  Therefore, 
Douglas  has  proposed  the  development  of  a  more  sophisticated  computerized  man 
BKidel  capable  of  being  Incorporated  into  comprehensive  programs  for  dose  cal¬ 
culations  within  realistic  vehicle  configurations  (i.e.,  SIQ<A).  Such  a  model 
would  provide  an  analytical  reference  against  which  dose  computations  based  on 
other,  simpler  models  could  be  measured.  Preliminary  analyses  have  Indicated 
the  development  of  a  Cartesian  coordinate  model  to  be  attractive.  In  this 
approach,  the  human  body  is  represented  as  a  number  of  separate  regions,  i.e., 
head  and  neck,  torso,  upper  arm,  forearm,  etc.,  each  region  be?ng  characterized 


by  a  rectangular  solid  geometry.  Superimposed  on  each  region  is  a  Cartesian 
coordinate  grid  structure  defining  an  array  of  grid  elements,  each  with  a 
specified  average  material  composition.  Advantages  of  such  a  model  relative 
to  extant  man  models  Include: 

1)  A  geometric  framework  for  which  grid  element  composition  data  may 
easily  be  prepared  from  anatomical  cross  section  drawings. 

2)  An  acceptable  approximation  to  the  spatial  distribution  of  body 
tissues  for  any  body  position,  by  rotation  and/or  translation  of 
the  several  separate  regions  with  respect  to  one  another. 

3)  A  geoBtttrlc  framework  providing  for  efficient  computation  of 
material  distributions  and  path  length  segments  along  rays  tra¬ 
versing  the  grid  elements  of  body  regions. 

4.2.2  Shield  Mass  Spatial  Distribution 

Douglas  has  Investigated  procedures  for  optimization  of  shield  systesm  associated 
with  nuclear  reactor  sources  on  board  space  vehicles  over  the  past  two  years. 

This  work  has  encompassed  the  developsient  of  mathematical  optimization  techni¬ 
ques,  the  Identification  of  shield  system  geometric  variables,  the  specification 
of  multiple  radiation  level  design  criteria,  and  the  automatic  computation  of 
radiation  levels  and  the  derivatives  of  radiation  levels  with  respect  to  the 
specified  variables.  Ihese  investigations  have  culminated  in  a  series  of  shield 
syst  m  optimization  programs  designated  "PATCH"  after  the  optimization  procedure 
employed.  Ihe  PATCH  technique  has  used  the  characteristic  thickness  of  designated 
shield  regions  as  variables  and  has  used  the  various  radiation  level  design 
criteria  to  constrain  the  values  of  these  thicknesses  during  the  optimization 
operations.  The  inawhematlcal  optimisation  is  based  on  the  formulation  of  the 
constraining  relationships  in  a  particular  manner  which  promotes  am  efficient 
iterative  solution  of  the  Lagrange  multiplier  representation  of  the  design  problem. 


Ihe  extension  of  these  optlioization  investigations  to  include  space  radiation 
sources  and  crew  conpertment  shielding  geometric  variables  is  now  in  progress. 

The  current  work  is  directed  toward  the  identification  of  a  flexible  framework 
of  variables  and  constraints  such  that  a  great  variety  of  space  vehicle  configur¬ 
ations  can  be  treated.  Therefore,  the  Initial  emphasis  is  on  use  of  parametric 
dose  data  from  SIGMA  computations  to  establish  functional  relationships  between 
dose  values  and  shield  thickness  variables,  as  required  for  optimization  opera¬ 
tions.  Completion  of  this  Initial  phase  of  the  shield  optimization  investigations 
will  serve  to  define  the  essential  feat  ^es  of  a  single  comprehensive  computer 
program  for  performing  dose  computations  in  generalized  quadric  geometry  and  for 
optimizing  the  shield  mass  distribution  subject  to  the  dose  design  criteria.. 

Specific  consideration  is  now  being  given  to  influencing  factors  such  as: 

1)  Crew  compartment  internal  arrangement,  especially  the  location  of 
bunks,  consoles,  rest  areas,  biowell,  etc. 

2)  Placement  of  equipment  and  stores  within  the  crew  compartment. 

3)  Division  of  shield  material  among  such  potential  locations  as  crew 
compartment  wall  structure,  blowell,  and  localized  and  personeQ 
portable  shields. 

*♦)  Shield  shaping  over  extended  shield  material  areas,  such  as  wall 
structure . 

5)  Time  dependence  of  dose  receptor  locations,  as  influenced  by  crew 
work-rest  cycles. 

An  additional  aspect  of  this  work  is  a  preliminary  treatment  of  uncertainties 
in  the  shielding  analysis.  This  phase  is  limited  to  the  translation  of  assumed 
uncertal!  les  in  the  intensity  of  the  incident  space  radiation  sources  into 
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uncertainties  either  in  the  dose  for  a  fixed  shield  design,  or  in  the  shield 
design  (particularly  vei^t)  for  specific  dose  design  criteria. 

U,3  Results  froB  SI(a<A  CoBputations 

Cosqputed  dose  data  for  varioiis  manned  orbiting  space  station  configurations  have 
been  prepared  using  the  OGRE,  CHARGE,  and  SIOtA  computer  programs.  *1110  depend¬ 
ence  of  dose  values  on  such  parameters  as  shield  thickness,  orbit  altitude,  orbit 
inclination,  and  shield  material  has  been  published  in  References  l6  and  17. 

An  example  of  these  dose  data  is  presented  here  for  the  typical  orbiting  space 
vehicle  depicted  in  Figure  9*  this  analysis,  doses  to  the  lens  of  the  eye, 
to  the  skin,  and  to  the  blood-forming  organs  vere  computed  for  an  astronaut 
located  in  the  center  of  the  laboratory  area.  The  man  model  used  in  the  calcu¬ 
lation  vas  comprised  of  two  elliptical  cylinders,  as  discussed  in  Section  4.2.1. 
Typical  results  for  dose  to  the  lens  of  the  eye  are  iUustwited  in  Figure  10. 
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